Dephasing at Low Temperatures 
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We discuss the significance and the calculation of dephasing at low temperatures. The particle 
is moving diffusively due to a static disorder configuration, while the interference between classical 
paths is suppressed due to the interaction with a dynamical environment. At high temperatures 
we may use the 'white noise approximation' (WNA), while at low temperatures we distinguish the 
contribution of 'zero point fluctuations' (ZPF) from the 'thermal noise contribution' (TNC). We 
study the limitations of the above semiclassical approach and suggest the required modifications. In 
particular we find that the ZPF contribution becomes irrelevant for thermal motion. 



The application of semiclassical considerations into the 
analysis of interference and dephasing is as old as the his- 
tory of quantum mechanics. A particular interest is to 
apply these considerations to the theory of diffusing elec- 
trons in a metal Wa. Dephasing, or the loss of 'phase' 
information, is the consequence of the interaction with 
some other environmental degrees of freedom, or with 
some measurement device. The widely accepted dogma 
is that dephasing is associated with leaving a trace in the 
environment. Lately this dogma has been challenged ex- 
perimentally || as well as theoretically Q, leading to 
an intensive debate associated with the question 

whether ZPF of environmental modes may lead to de- 
coherence at the limit of zero temperature. The most 
physically-appealing theoretical considerations are based 
on the application of Feynman- Vernon (FV) formalism, 
and simple semiclassical considerations. In particular, 
the zero temperature decoherence found within the ex- 
actly solvable Caldeira-Leggett (CL) model, is most puz- 
zling. Thus, it is of fundamental importance to address 
the following: (1) What are the limitations of the semi- 
classical strategy; (2) What is the applicability of results 
that are based on the CL model; (3) What are the limita- 
tions of the 'one particle' picture when applied to a many- 
body system. In this letter we are going to explore these 
questions systematically, using a well controlled strategy. 
Under consideration is the motion of a particle under the 
combined influence of static disorder and a dynamical en- 
vironment. Our reasoning consist of the following three 
stages: (a) We take the static disorder into account and 
write the transport probability amplitude as a discrete 
sum over classical trajectories; (b) We take the stochas- 
tic nature of the environment into account by consider- 
ing the influence of an effective stochastic potential; (c) 
We take the full dynamical nature of the environment 
into account by using the FV formalism. Having an en- 
vironment whose temperature is T, and a particle that 
is injected with an energy E (not necessarily thermal), 
one should make a distinction between the cases of high 
and low temperature T, and analogous distinction be- 
tween the cases of large and small energy E. We shall 
argue that the validity of the semiclassical approach is 
not restricted to high bath-temperatures. However, at 



low bath-temperatures, an essential modification is re- 
quired in case of motion with small energy E. 

It is assumed that the motion of the particle under 
the influence of the static disorder is diffusive. The 
transport probability amplitude can be written as a sum 
^2 A a exp(iS[x a ]/fi) over classical trajectories that con- 
nect the observation point with the injection point. We 
shall denote by t the total time of the motion, and S[x a ] 
are the corresponding action integrals. The statistical 
properties of the differences Xbfa) — x a (ti), where a&6 is 
a pair of trajectories, will play a major role in later cal- 
culations. These statistical properties can be taken into 
account via a single function P(k,u>) that reflects the 
power-spectrum of the motion j7|. For diffusive motion 
a practical approximation is 



P(k,u) = 2Dk 2 I {{Dk 2 + 1/tf + u 2 ) 



(1) 



where D is the diffusion coefficient. In the large k regime 
where k>(Dt)~ 1 / 2 this function is a properly normalized 
Lorentzian. Note that the collision frequency lu~v 2 /D, 
where v is the velocity of the particle, should be used as a 
cutoff to the slow 1 juj 2 power-law decay. Beyond k^v /D 
the above expression is not valid and the power spectrum 
of the motion is ballistic-like. 

We now take the stochastic nature of the environ- 
ment into account by introducing into the Hamiltonian 
a stochastic potential that satisfies 



(U(x" ,t")U(x' ,0) = <t>{t"-?) ■ w(x"-x') 



(2) 



It is assumed that w(r) depends only on \r\. The inten- 
sity of the noise is characterized by the parameter 



4>{t)oIt ■ \w"(0)\ 



(3) 



The power spectrum of the noise </>(w) is the Fourier 
transform of <P(t). We shall assume ohmic environment, 
meaning that at the classical limit 4>{oj) = v up to some 
cutoff frequency 1 /r c which is assumed to be larger than 
any other relevant frequency scale. Thus, in the classical 
limit we can use the WNA, namely 4>(t) = v8(t). For 
the quantum mechanical see ([lO]) . Without loss of gener- 
ality we shall assume the normalization w"(0)— — 1. The 
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d-dimensional Fourier transform of w(r) will be denoted 
by w(k). The mode-density (after angular integration) is 
g(k) = {C d /{2Tr) d )k d ~ 1 w{k) 1 where C d is the total solid 
angle. We shall assume that 



j(k) = C£ 2 



for k < 1/ 



(4) 



where £ characterize the spatial scale of the correlations, 
and C is a dimensionless constant. In case of short range 
Gaussian-type correlations a equals simply d. For the 
long range Coulomb interaction to be discussed later it 
equals d— 2. In order to have a well defined model we 
must have |u/'(0)| < oo therefore only —2 < a is mean- 
ingful. The regime — 2<cr<0 is well defined but it requires 
special treatment since w(0) diverges. 

The path-integral expression for the probability to 
propagate from the injection point to the observation 
point constitutes a double sum J J Vx'Vx" over the path 
variables x'(t) and x"(t). Averaging over realizations of 
U one obtains the influence functional exp(— Sn /h 2 ). See 
H for details. The noise action functional is 

S N [x',x"} = ~ / [ dtidh <j>(te-ti) x 
^ Jo Jo 

[w(x%-x") + w{x' 2 -x' x ) - 2w(x'^-x[)} (5) 

where shorthand notation for x(ti). The corre- 

sponding semiclassical expression for that probability is 



a b 



A a A* b exp 



S N [x a , x b ] 



.s[x a }-s[x b y, 

exp [ i ) (6) 



It is obvious that the interference contribution (the terms 
a^b) is suppressed due to the noise, while the classical 
(diagonal) contribution survives f. This is the the de- 
phasing effect. See general discussion in We are 
interested here in the computation of these dephasing 
factors, as well as in illuminating their physical signifi- 
cance. For the purpose of calculating the typical value of 
the dephasing factor, (H) can be replaced by 



(S N ) « t / g{k)dk 

JO JQ 



dio 



(Lu)P(k,L0) 



(7) 



If the integral on the right hand side is independent of t, 
then the typical dephasing factor in (^) can be written 
as a simple exponential exp(— t/r<p). 

The full analysis should take into account the dynam- 
ical nature of the environment. Namely, one should con- 
sider an Hamiltonian H = Ho(x,p) + H env (x, Q a , P a ), 
where the latter term incorporates the interaction with 
environmental degrees of freedom. It is possible in princi- 
ple (but generally not in practice) to use the FV formal- 
ism in order to derive an appropriate influence functional 
exp(z5p/7i — Spi/h 2 )- The Fluctuation-Dissipation theo- 
rem (FDT) implies that if Sn is known, and the tempera- 
ture of the bath is further specified, then also some of the 



dissipative properties of the environment are determined 
uniquely. Therefore it is plausible that the CL procedure 
of constructing an effective harmonic-bath, is useful in 
order to derive an actual expression for the friction func- 
tional Sp. Indeed, this strategy has been adapted in ||] 
and lead to the introduction of the DLD model \. The 
interaction with the bath-oscillators is 



Hi = ^c a Q a u(x-x a ) 



(8) 



Here x a is the (fixed) location of the a oscillator and Q a 
is its dynamical coordinate. The bath-oscillators are dis- 
tributed uniformly all over space. The interaction of the 
particle with each of the oscillators is described by u(r). 
The range of the interaction is £, and c a are coupling 
constants. For an ohmic bath the following expression 
(generalized here for any dimension) has been derived: 



S F = 



r\ I dr Vw(r) ■ R 
o 



(9) 



where rj is a friction parameter, and the path variables 
are r = x"-x' and R = {x"+x')/2. From FDT it fol- 
lows that if an ohmic environment is characterized by 
a temperature T then the friction parameter should be 
i] = ^/(2fceT) We shall assume an environment that is 
characterized by a short spatial autocorrelation scale £, 
such that the classical trajectories are well separated with 
respect to this microscopic scale. Under such circum- 
stances it has been observed in || that Sf will have no 
effect on the interference contribution. This statement 
does not hold in case of the CL model. The CL version for 
Sf is obtained by taking in ([)]) the limit £ — > oo, which is 
equivalent to the formal substitution w(r) = —r 2 /2. Av- 
eraging the factor exp(iSp/h) over diffusive trajectories 
one obtains, as in ||, a non-generic factor exp(— (t/T v ) 2 ) 
where \Jr lp = r/D/h. This particular result turns out to 
be identical, up to a logarithmic factor, with the gen- 
uine result (|l2|), to be discussed later. However, it is not 
consistent to use the CL version for Sf in the present cir- 
cumstances, and therefore the approach of || does not 
apply. With the above observations, our semiclassical 
strategy implies that Sf of a generic environment has no 
consequence on the analysis of dephasing, and Eq.(^) is 
still valid. Our main conclusion below will be that this 
(semiclassically-based) statement fails for a low energy 
particle. 

We turn now to discuss some actual results for the de- 
phasing rate. For ohmic bath the symmetrized power 
spectrum of the noise is 



<P(uj) = T)\cj\ h coth 



h\u)\ 
2k B T 



for |w| < l/r c (10) 



For high temperatures l/r c < ksT/h, one can use the 
WNA. Substituting ^{r)=2r}k B T into (§), one obtains 
H the universal high temperature result 



2 



2nk B T£ 2 



for < a 



(11) 



This result does not apply for — 2<er<0 (for elec- 
trons d<2), because w(0) diverges. Still, using the 
WNA and doing some simple manipulations, one obtains 
S N =2nk B T£ 2+a (Dt)- a / 2 t, leading to a dephasing factor 
of the type exp(-(t/r (p )( 2 -' T )/ 2 ), where 1/t^_oc T 2 ^ 2 -^ 
in agreement with the well known results 




FIG. 1. The (k,u>) plane. The shaded regions indicate those envi- 
ronmental modes that are effective in the dephasing process. The 
darker region indicates a possible excess contribution due to ZPF. 

We are interested in going beyond the WNA. It is use- 
ful to use (Q) in order to perform the actual calculation. 
The domain of integration is illustrated in Fig.l. If the 
temperature is high enough, such that D/£ 2 < k B T/U, 
then one can still use the WNA result. As the temper- 
ature becomes low, such that k B T/h < D/£ 2 , the de- 
phasing rate 1/t v becomes larger than the value which 
is predicted by the WNA. This is due to the ZPF in the 
frequency zone k B T/h < lu where <j>(uS)=TiriU) . Later we 
shall see that E/h should be used as a cutoff for the u> 
integration. Therefore the present analysis, and ( fH[ ) be- 
low, does not apply to a low- energy particle. Considering 
the ZPF contribution, the integral in (Q) is dominated 
by (k,uj) modes that are concentrated along the curve 
uj = Dk 2 . Moreover, most of the contribution comes from 
modes with large wavenumber, namely fc~l/£ Thus the 
ZPF contribution to the dephasing rate is essentially as 
in H] and applies at low temperatures to any — 2 < a 



C 



V 3 / ZPF 



(2+o> 



In 1 




x -nD . (12) 



At low but finite temperatures, meaning l/t < k B T/h <C 
D/£ 2 , one should consider the TNC that comes from the 
lower shaded region of Fig.l. For — 2<er<0 the TNC de- 
phasing factor is determined by k^(Dt)^ 1 ^ 2 modes. For 
these modes Dk 2 <k B T ' j%, and therefore the previously 
discussed WNA applies. On the other hand for 0<er one 
obtains 



C'£ 2+ ° 



k B T 
TiD 



a/2 



2 V k B T 



(13) 



where C is a numerical factor. The above expression is 
valid for 0<a<2, where the dephasing process is dom- 
inated by modes with k~(k B T /UD) 1 / 2 . For 2<a the 
dephasing process is dominated by modes with k^l/£, 
and (|l^) should be modified by the replacement a — s- 2. 

The FV path-integral expression is exact in principle. 
On the other hand, our semiclassical expression (JsJ) in- 
volves the stationary-phase approximation. Therefore it 
is important to understand, physically as well as math- 
ematically, the validity limits of the semiclassical strat- 
egy. Each stationary-phase point of the exact FV path 
integral is a pair a&cb of real classical trajectories that 
correspond to the motion of a fictitious classical test par- 
ticle. Let us assume for simplicity that akb are loops 
related by time-reversal. It is also essential in the follow- 
ing argumentation that the interaction with the environ- 
ment is short range, as in the 0<cr DLD model with £ 
which is a small scale. It follows from the definition of 
the influence functional that under such circumstances 
P = 1— cxp(— Sn[o]) is the probability for a fictitious test 
particle to leave a trace along the way (i.e. to change the 
quantum mechanical state of at least one bath-oscillator 
along the loop). For simplicity one may consider the 
restricted problem of a particle that travels across a sin- 
gle bath-oscillator, meaning that (||) includes only one 
term. It is well known |)| that the semiclassical ap- 
proximation is equivalent to a self-consistent mean-field 
scheme, where it is assumed that the wave-function for 
the particle-oscillator system can be written as a prod- 
uct of a scattered particle and driven oscillator. Such 
ansatz implies that for weak scattering we can treat the 
particle as moving with constant velocity v and solve for 
the oscillator. It turns out that this reduction requires 
the assumption of small energy transfer. Therefore, one 
should anticipate problems once an oscillator with uj a 
larger than E is involved. In the latter case, there is no 
justification to think of the particle as decoupled from 
the bath, moving with some constant velocity, capable of 
exciting the oscillator along the way. Therefore the cor- 
responding probability P loses its physical significance. 
Still, the interference contribution of the specified loop 
will be suppressed. It is true that a low-energy particle 
(uj a < E) is not capable of exciting the oscillator, but 
there is a finite probability to be scattered elastically, 
and consequently the original interference contribution 
is suppressed. However, this suppression does not have 
the meaning of genuine decoherence. Rather, it reflects a 
coherent "re-normalization" of the bare dynamics. One 
wonders whether exp(— 5/v[a]) gives an estimation for the 
suppression of the original interference contribution due 
to these elastic scattering events. Unfortunately this is 
not the case. The elastic scattering off an oscillator in- 
volves a second-order process of virtual emission followed 
by absorption of a quanta hu a . Therefore the elastic 
scattering probability is proportional to c„ . At the same 
time P is proportional to c 2 in leading order. Therefore, 
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under such circumstances, exp(— Sjv[a]) does not reflect 
any physically meaningful quantity. It should be empha- 
sized that the above argumentation does not apply once 
£ becomes large compared with the distance that is ex- 
plored by the particle. In case of the CL model {£ — > oo) 
the decoherence process is no longer determined by the 
scattering mechanism. Rather, we have a crossover to 
a non-perturbative spreading mechanism. See detailed 
discussion of this point in |?J . 

The semiclassical strategy cannot be trusted if some 
of the effective bath-oscillators are such that E/Ti < 
uj a . For diffusive motion such a situation will occur if 
E/Ti < D/£ 2 . Note that the notion of large/small en- 
ergy is in complete analogy with the notion of high/low 
temperature. For thermal motion E ~ ksT and the 
two notions coincides. In order to see how expression 
(fjj) should be modified for small-energy motion, let us 
obtain it using an elementary perturbative calculation. 
This is straightforward for short range interaction (0<tr) 
since it is plausible that 1/t v has then the significance 
of inelastic scattering rate. The leading order inelastic 
scattering probability is expressed in terms of the diag- 
onal matrix elements of the scattering matrix, namely 
P = (1 — \ (nip\S\n^p)\ 2 )0. The subscript (3 implies ther- 
mal average over the quantum-mechanical states n of 
the bath-oscillators. The incoming particle is described 
by the wavefunction tp. Note that the precise definition 
of 'dephasing-rate' becomes vague once the semiclassical 
approach is abandoned. It seems plausible that P should 
be averaged over eigenstates of the disordered potential. 
Using second order perturbation theory 



(nip\S\nip) 



1 

1 o 



dt 2 dt 1 (mj}\Hi{t2)'rti{t 1 )\n^) 



o Jo 



Standard manipulations which are based on the defini- 
tion of the ohmic DLD model lead to the result P = 1/t v 
where 



1 



1 



dk 



(2?r) d J 2ir 



du; 



w{k)$(w) P qm {-k,-Lu) (14) 



Here = hrj\uj\n{uj) is the non-symmetrized version 

of </>(u>), where n(ui) = l/(exp(huj/kBT) — 1) for 0<ui, 
and n(u>) = (l+n(|w|)) for lo<0. At zero temperature 
<!>(0<w)=0. The quantum-mechanical power spectrum 
of the motion is defined as the (complex) Fourier trans- 
form of P qm (k, t) which is the correlator of the operator 
cxp(zA;a;). Using the semiclassical estimate P qm (k,oj) w 
P(k,u>) one obtains again the integral expression in (Q) 
for the dephasing rate. The validity of the semiclassi- 
cal estimate for P qm {k,uj) is restricted to the frequency 
range \u\ < E/%. Quantum- mechanics cannot support 
larger frequencies! For thermal motion the result of the 
quantal FDT is related to the result of the classical FDT 
as follows: 



P qm (k,u>) 



h\uj\ 



Thus, it is suggested that for thermal motion ksT/h 
should serve as an effective cutoff in the w integration 
of Eq.(0). Consequently the ZPF contribution for the 
dephasing rate should be omitted. For non-thermal mo- 
tion E/Ti is the proper cutoff and some (or all) of the 
ZPF contribution should be included. 

Finally, one wonders whether additional modifications 
are required once we turn to treat (semiclassically) the 
problem of dephasing of electrons in a metal, taking 
into account the presence of Fermi sea. Here we con- 
sider ballistic-like motion as a test case. The aver- 
age scattering rate can be calculated by using a ki- 
netic picture. The transition rate for unit volume 
is W(p'\p) — w(k)<fr(uj)/Ti 2 , where k—(p'—p)/h and 
oj=(E(p')—E(p))/h. The average scattering rate is pro- 
portional to / / dp'dp(l-f(E'))f(E)W(p'\p). Standard 
manipulation leads to Eq.(O) with (15). Note that for 



the ballistic-like motion under consideration one should 
substitute P(k, ui) = l/(i>|fc|) for |w|<i)|fc and zero other- 
wise. This calculation implies that the dephasing rate is 
similar to the inelastic scattering rate also for the many- 
body case (provided 0<a). 
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It is always true that < Sn- However, the concept of 
dephasing implies that Sn — ► oo for long trajectories. The 
assumption of having short r c is important here. If we take 
4>(t) = const, then some important interference contribu- 
tion will not be suppressed in the t — ► oo limit. See || 
for details. This observation is essential in order to un- 
derstand, within the path-integral approach, why static 
disorder leads to either weak or strong localization effect. 
DLD are the initials of "Diffusion Localization and Dissi- 
pation" . Obviously, these three effects comes out naturally 
from the solution of the DLD model. Note however that 
localization effect is out of the scope of the preset appli- 
cation. See the previous footnote in this respect. 
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